Nonlinear waves in a positive-negative coupled waveguide zigzag array 
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We consider the coupled electromagnetic waves propagating in a waveguide array, which consists 
of alternating waveguides of positive and negative refraction indexes. Due to zigzag configuration 
there are interactions between both nearest and next nearest neighboring waveguides exist. It is 
shown that there is a stop band in the spectrum of linear waves. The system of evolution equations 
for coupled waves has the steady state solution describing the electromagnetic pulse running in the 
array. Numerical simulation demonstrates robustness of these solitary waves. 



PACS numbers: 42.65.Tg, 42.79.Gn , 42.81. Qb 



I. INTRODUCTION 

Recently new applied physics field is referred to as trans- 
formation optics develops [T-'B^. The transformation op- 
tics owes its origin to the artificial materials — metama- 
terials, which are featured by unusual electrodynamics 
properties [THIQ]. For instance these media are charac- 
terized by a negative refraction index (NRI) when the 
real parts of the dielectric permittivity and the magnetic 
permeability are simultaneously negative in a certain fre- 
quency range. 

The negative refraction in metamaterials manifest 
themselves when the wave passes through the interface 
between such medium and a conventional dielectric or 
in the case of the surface waves [11] [12]. The novel 
kind of nonlinear interaction of the guided waves can 
be realized in a nonlinear oppositely-directional coupler 
\T3\ fT4] . This coupler consist of two tightly spaced non- 
linear/linear waveguides. The sign of the index of re- 
fraction of one of these waveguides is positive and the 
index of refraction of other waveguide is negative. The 
opposite directionality of the phase velocity and the en- 
ergy flow in the NRI waveguide facilitates an effective 
feedback mechanism that leads to optical bistability [T5] 
and gap soliton formation [TH [1^ . 

The optical waveguide array provide a convenient 
setup for experimental investigation of periodic nonlinear 
systems in one dimension [17 . Nonlinear optical waveg- 
uide arrays (NOWA) are a natural generalization of non- 
linear couplers. NOWA with a positive refractive index 
have many useful applications and are well studied in 
the literature (see for example [T8H2Q] ). If the sign of 
the index of refraction of one of waveguides in NOWA is 
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FIG. 1: A schematic illustration of zigzag coupled waveguides 
array (a); and zigzag positive-negative coupled waveguides array 
(b). 



positive and the index of refraction of other neighboring 
waveguide waveguide is negative the alternated NOWA 
will be obtained [23'-'25 . 

Usually the coupling between nearest neighboring 
waveguides is taken into account. It is correct approxima- 
tion for strong localized electromagnetic wave in waveg- 
uide. However, the coupling between both nearest neigh- 
boring waveguides and the next nearest neighboring ones 
can be introduced by the use of a zigzag arrangement 
[2T| I22] . Let 1^1) is an angle between the lines connect- 
ing the centers of neighboring waveguides. In a linear 
array this angle is tt. In a zigzag like array at '^5 ^ 7r/2 
the coupling between the nearest neighboring waveguides 
and the next nearest neighboring ones is approximately 
the same. Nonlinear optical waveguide zigzag arrays 
(NOWZA) can be considered as generation of NOWA. 

In this paper we will consider nonlinear optical waveg- 
uide array that looks like zigzag array discussed in [21]. 
However in the case under consideration the refractive in- 
dexes of the neighboring waveguides have a different sign 



(Fig. 1 b). The dispersion relation of the hnear wave 
win be determined. We suppose that only waveguides 
with PRI are manufactured from cubic nonlinear dielec- 
tric. Steady state solitary waves of the particular kind 
will be found. Numerical simulation of the steady state 
waves collision will be used to illustrate the robustness 
of these solitary waves. 



II. THE MODEL FORMULATION 

Let us assume that waveguide marked by index n is char- 
acterized by positive refractive index (PRI), the nearest 
neighboring waveguides with indexes n — 1 and n + 1 pos- 
sess negative refractive index. If the electromagnetic ra- 
diation is localized in each waveguide the coupled wave 
theory can be used. In the case of the array which is 
prepared in the form of line (i.e., all centers of waveg- 
uide placed on one line), only the interaction between the 
nearest neighboring waveguides is essential. In this con- 
figuration the angles between the lines connecting waveg- 
uides are equal tt. However, if the array is deformed in 
the form of zigzag, where the angles between the lines 
connecting waveguides are equal approximately 27r/3, in- 
teraction between both the nearest neighboring and the 
next nearest neighboring waveguides will be important 

m- 

The system of equations describing the wave propa- 
gation in a nonlinear zigzag waveguide array with alter- 
nated sign of refractive index reads as 

M ^ + ^ I en + ci(en+i + en-i) + 



+C2(en+2 + en-2) + ri\enfen = 0, 

d d\ 

1 I -^ + ^ I en+1 + ci(en+2 + en) + 
+C3(en+3 + en-i)r2|enpen = 0, 



(1) 



(2) 



where e^, is the normalized envelope of the wave localized 
in n-th waveguide. Coupling between neighboring PRI 
and NRI waveguides is defined by parameter ci . The pa- 
rameter C2 (C3) is coupling constant between neighbor- 
ing PRI (NRI) waveguides. The all functions en{C->^)i 
independent variables C? '^ ^^^ other parameters are 
expressed in terms of the physical values represented in 
[23 . We will assume that NRI waveguides are linear 
ones (r2 = 0). The configuration of these alternating 
waveguides will be remarked as asymmetric alternating 
nonlinear optical waveguide zigzag arrays (ANOWZA). 



III. LINEAR PROPERTIES OF THE 
ALTERNATING WAVEGUIDE ZIGZAG ARRAY 

let us consider the case of the linear waveguides combined 
into zigzag alternating waveguide array. This optical sys- 



tem is described by following equations 
d d 



dC dr 



Cn +Ci(en+1 +en-l) + 



+C2(en+2 +en-2) = 0, 
d d\ 

^ ^ ~Bt ~ 8~ I ^^^^ ~ ei(en+2 + en) 

-C3(en+3 +en-l) = 0. 



(3) 



(4) 



To find the linear wave spectrum we can employ the 
presentation of the envelopes in the form of harmonic 
waves 

Substitution of this expression in ^ and Q leads to 
a system of the algebraic linear equations respecting A 
and B. Nonzero solutions of these equations exist if the 
determinant 



det 



(7 -I- cj -I- 73 71 

-71 (7 - cj - 72 



is equal to zero. Here the parameters 71 = 2ci cos (/), 
72 = 2c2 cos 2^, and 73 = 2c3 cos 2^ were introduced. 
This condition results in equation 



where 



(cj + cjo)^ =7i + {Q-Qof 



2ujo = (72 + 73), 2(70 = (72 - 73). 



Thus, we obtain the dispersion law for linear waves in 
linearized ANOWZA 



00 



{q) = -ujo ± y 7i + (<7 - ^0)^ 



i^)(n\ = 



(5) 



Expression ([5| shows that (a) the forbidden zone (gap) 
in spectrum of the linear waves exist Auj = 271 , (b) spec- 
trum is shifted along both frequency axis and wave num- 
bers axis, (c) the form of spectrum likes the spectrum 
for linear oppositely-directional coupler jl4| [23l [25] . The 
gapless spectrum appears only when condition (j) = 7r/2 is 
hold. In this case the radiation propagates along waveg- 
uides with the same refractive indexes. Energy fiux be- 
tween neighboring waveguides is zero. 



IV. NONLINEAR WAVES IN ANOWZA 

The simplest approximation is based on following ansatz 

e„(C,r) = A(C,r)e^*", e„+i(C,T) = i?(C,r)e^*("+i), 

where A and B are the envelopes of the quasi-harmonic 
waves. It allows to do reduction of the equations ([l}-([2]) 
and to obtain the equations 



d_ d_ 

d__d_ 
dC dr 



A + -iiB + -(2A + r\A\^A = 0, (6) 
B--iiA- 73B = 0. (7) 



At ^ = 7r/2 this system of equations is splitting in 
two independent equations. At ^ = 7r/4 ([gl-Q is trans- 
formed to that, which was considered in pSl fM]. 

Real amplitudes and phases can be introduced by 
means of the formulae A = aiexp(z(^i) and B = 
a2 ex.-p{i(p2) . The system of equations (|6| and ^ can 
be represented in the real form: 
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dr 



ai 


= 7ia2sin$, 


a2 


= 7iai sin<l>, 


^1 


<^2 

= 7i — cos <l> + 72 + 


^2 


= -71^ cos ^ -73, 

^2 



where ^ = cpi — ip2. 



A. Steady state w^aves 

Among the different kinds of the waves the steady state 
waves attract attention. These waves correspond to solu- 
tion of the wave equation depending only one particular 
variable. Let this variable be 



71 






where /3 is parameter that defines the group velocity. If 
we introduce new envelopes ui and U2^ according to for- 
mulae Ui = yi + /^<^i '^2 = VI ~ /^<^2, then -ui, U2 and ^ 
are governed by following equations 



dui 

dU2 



= ^i2sin$, 



= -Uisin^, 



J., / ^1 , ^2 \_^ , ^,^.2 



U2 Ui 



COS^ + ^iil, 



(8) 

(9) 

(10) 



where 




From equations ([8| and ([9| the first integral of motion 
follows 



U2 = Ci = const. 



(11) 



If both parts of (ITO ) multiply on U1U2 cos ^ and take into 
account (|8| andTO), then the second integral of motion 
can be obtained 



U1U2 cos$ 



i9 



C2 



const. 



(12) 



In depending on the boundary conditions the different 
solutions of the equations ([8|-( 10 ) can be found. Any of 
these solutions will be describing the steady state waves. 



B. Solitary Avaves 



Solutions of the equations M-(IO) under boundary con- 



ditions iii^2 ^ at 1^1 ^ 00 correspond with steady state 
solitary waves. According to the boundary conditions the 
constants Ci and C2 of (11) and (12) are zero. Hence 



ui -ul=0, £ cos <l> + - + -ii? =0, e = ±1. 

I z 1 2 4 



The system of equations (p|)-(10) with taking into ac- 



count U2 — eu\ can be reduced to the following one 
du\ 



eu\ sin<l>, 



Taking into account the second integral of motion the 
equation for i^i is represented as 



du\ 



= ul 



'& 



Substitution ui = w ^/^ and some little manipulations 
leads this equation to 



-^ ) ='iA^{w-Wi){w + W2), 



where 
A2= (1 



4 



Wl 



^ 



4(1 -V2)' 



W2 



^ 



4(1 + V2)* 



Integration of the last equation for w can be done (see 
[14 ). It results in 

w = wo^ws cosh[2A(^ - ^0)], 

where ^0 is a constant of integrating. Thus, we have 

61^ 1^1 



Wo 



4(1-^74)' 



Ws 



4(1 -(^74)* 



In term of original variables the functions 1^1^2(0 can 
be written as 



uliO = uliO 



4AV|^| 



cosh[2A(e-eo)]+V2' 



(13) 



So, the real envelopes of the solitary wave propagating 
in ANOWZA are presented by following expressions 



alio 



4A2 



-, (14) 



|t9|(l+/?){cosh[2A(C-Co)]+5/2}^ 

2 4^^ 

'''^^^ " M(l-/3){cosh[2A(e-eo)] + V2}' ^^^^ 



The phase difference i> evolves according to expression 

$(0 = $(-^) + sgn(^)§(V2; Xo), 

where auxiliary function 

(1 -a^)^/^e^o 

S(a; Xo) = arctan ^z , 

1 + ae^o 

and the argument Xq = 2A(^ — ^q) were introduced. 

Due to that at ^ ^ — oo the derivative dui/d£, to be 
positive, the value of the <!>(— oo) has to equal 67r/2, where 

£ = ±1. 

C. Algebraic solitary waves 

The solutions of the equations ([8|-([l0| describe the ex- 
ponentially decaying wave fronts. However some times 
the solitary waves can be decreasing as ^ 1/^^. 

The solutions ( 14 )-( 15 ) are characterized by parameter 
A, which can be set to zero in limit of \5\ -^ 2. Since, 
the solutions found here are correct if —2 < S < 2. How- 
ever, on the boundaries of this interval we have to refine 
behavior of the solitary waves. Let us consider the func- 
tion u^{^) near little values of A. At A <C 1 from (13) it 
follows 



_ 4 
_ 4 
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4AVt? 



■2A2(e- 
(1 



■ eo)2 + 
-<5V4) 



5/2 



■5/2- 



■2(1 
(1- 



6/2){l + 5/2) 



(1 



■6/2) 
1 



-2(1- 

S/2 



.5/2)il + S/2)i^-^oY 



l + 2{l-5/2){^-^oy 



Near left boundary of the interval for permissible value 
of parameter J, the function i^^(^) can be found using the 
substitution S = —2 -\- e^ where £ <C 1. It is follows 



\o 



hm - 



2-£/2 



l + 2(2-e/2)(^-Co)^ 



Near right boundary of the interval for permissible 
value of parameter S by the use the substitution 5 = 2 — e. 
It results in following expression 

u2(e) = hmlf— -^^?-— ^')=0. 

e^O § \1 + S{^ - ^o)y 

Thus, when parameter 5 runs to —2 the solutions of the 
equations ([8l)-(|10|) take the form of the algebraic soliton 



am = 

4i0 = 



^(l + /3i){l + 4(e-eoP}' 



(16) 
(17) 



t?(l-/?i){l+4(C-6)2}- 

Here /3i corresponds with S = —2. If 5 runs to +2, the 
amplitudes of the waves |8|-(10) are equal to zero. 



V. LIMITATION OF THE VELOCITY 
PARAMETER /3 

The velocity parameter (3 is limited due to condition \S\ < 
2. With the definition 




1-/3 
we can introduce function 

71 
depending on 



72^ 

71 m' 




fi 



It is odd function of ja. Hence the condition \S\ < 2 
will be held in intervals ji- < ji < fij^ and — /i_ > (jl > 
— /i+, where /i_ and /i+ are positive roots for equation 
|F(/i)| = 2. For positive /i if minF(/i) < 2 the real roots 
exist. We can find that the function F^ji) attains ex- 
tremum at /i = ±\/72/73- If /^ is positive value then 
this extremum is minimum, else if ji is negative one the 
extremum is maximum. The value of the function F^fi) 
dX ji = +V72/73 is 2y^7273. Thus the solitary waves are 
exist if y^7273 < |7i|. From the definitions of parame- 
ters 71,2,3 it results in 7273 = 4c2C3 cos^ 2(j). In discussed 
model the coupling coefficients ci,2,3 can be made posi- 
tive ones. Hence, the product 7273 is positive too, or is 
equal to zero at = 7r/4. 

Parameter 5 can be represented as 




It follows that the phase shift for complex envelopes in 
neighboring waveguides defines the sgn(5. For example 
parameter 5 will be negative one in intervals tt/2 < (p < tt 
and 37r/2 < ^ < 27r. 

If /i > the roots of equation F(/i) = 2 take the fol- 
lowing form 



Z^- 



M+ 



N 
I72I 

N 
I72I 



1 



1 



7273 ^ 
7273 ^ 

7? ; 



The boundary of interval for permissible values of the 
velocity parameters /3± are determined by the expression 



P± = 



f^l 



1 



M^ 



1 



For example, let us assume that 72 
case we have 



73 = 0,571. In this 



/i+ ^ 3, 73 /3+ ^ 0, 866, 
/i_ ^0,266 /3_ ^-0,865 





FIG. 2: Collision of two steady state pulses in the case of 72 
73 = 0. 





FIG. 3: Collision of two steady state pulses in the case of 72 
73 = 0.01. 



It should be pointed that if the next nearest neighbor- 
ing waveguides are not coupled as in [14], then steady 
state waves exist for \B\ < 1. 



VI. NUMERICAL SIMULATION OF THE 
SOLITARY WAVES COLLISIONS 

To demonstrate stability of the solitary waves corre- 
sponding with the system of equations ([6| and (JTl) the 
collision between two solitary waves has been simulated. 
The real envelopes were taken as (14 andp^Sl). The real 



phases corresponding them were chosen as following ones 






sgn(i9)S(0,2O), 

; + isgn(^)S(0,20), (18) 



where Q = 71 [Cj — /3(r — rj)]{l — /3^)-^/^ The arguments 
Q Tj define the initial position of the j-th solitary wave. 
To produce a collision between two solitary waves we 
used solutions of the equations ^ and (jl with param- 
eters P = 0.4 and P = —0.4 as initial conditions for 
these equations. The pulse with f3 = +0.4 was located 
at C = 50 and pulse with /3 = —0.4 was located at C = 0. 
The coupling constant 71 was set as unite. For simplicity 
the coupling constants 72 and 73 are assumed equal, i.e., 

72 =73. 

The case of the linear chain of the alternating refrac- 
tive indexes waveguide array corresponds to waveguide 
system with (72 = 73 = 0). Collision between two soli- 
tary waves is represented in Fig J2] (see also [14 ). Figures 
show that collision is elastic. 

We found that the collision between counter propagat- 
ing pulses is elastic for the coupling constants 72 = 73 
that is taken from interval [0.001, 0.0075]. Little radia- 
tion appears where the coupling constants are taken from 
interval [0.0075, 0.01] Figjs] For the coupling constants 
are from interval [0.001, 0.021] the amplitudes of the ra- 
diation are increasing. 

Increasing the coupling constants 72 = 73, that are 
belong to interval [0.0215, 0.03], leads to the transferred 
wave decreasing as it is shown in Fig. |4] 

When 72 = 73 > 0.02 the reflected wave appears in 
NRI waveguide as a result of the reflection of the incident 
solitary wave (/3 = —0.4) from solitary wave (/3 = 0.4) 
propagated in the opposite direction in ANOWZA. The 





FIG. 4: Collision of two steady state pulses in the case of 72 
73 = 0.022. 



reflected wave amplitude grows up to some value and 
thereafter it decrease in interval [0.025, 0.05], Figjsland 
FigJ6] The reflected wave is disappeared at 72 = 73 ^ 
0.075. 

If the coupling constants belong to interval 
[0.08, 0.135] the steady state solitary waves are 
akin to elastic interacting waves. There is no radiation 
after collision, however the velocities of resulting waves 
vary considerably. See Fig. [7|and Fig. [2j 

If 72 and 73 are taken from interval [0.138, 0.15], the 
weak linear wave in NRI waveguide is generated. This 
linear wave is absorbed by nonlinear transferred wave 
when 72 and 73 rank among the interval [0.151, 0.16]. 
Reflected wave appear again when 72 and 73 are among 
the intervals [0.17, 0.2] and [0.31, 0.34], see Fig. [s] 

The steady state pulses collision does not generate the 
radiation if the coupling constants 72 = 73 are taken 
from interval [0.21, 0.3]. However the velocities of the 
scattered pulses can be strongly varied with respect of 
initial values. 

Also we had considered the interaction of the inci- 
dent solitary wave corresponding to /3 = 0.4 with quasi- 
harmonic wave. Initial incident pulse was located at 
C = 50, the quasi-harmonic wave was generated at C = 0. 
This wave is described by following expression 

/(r) = /o[tanh(r - 40) - tanh(r - 70)] sinuJbgT. (19) 





FIG. 5: Collision of two steady state pulses in the case of 72 
73 = 0.06. 







FIG. 6: Collision of two steady state pulses in the case of 72 
73 = 0.08. 



FIG. 9: Steady state pulses crossing the quasi-harmonic plateau 
like wave, 72 =73 = 0. 





FIG. 7: Collision of two steady state pulses in the case of 72 
73 = 0.13. 



where Ui^g is the frequency of modulation, /o is the am- 
phtude of this plateau like pulse. Numerical simulations 
have been proposed for 0.1 < /o < 3 and 0.7 < uj^g < 3. 
It was found that initial steady state solitary wave is 
extremely prone to damage if the coupling constants 72 = 
73 are less than 0.05. FigJ9] shows the propagation of 
incident pulse trough quasi-harmonic wave (19) at /o = 
1.5 diiid ujbg 



0.7 and at 72 =73 =0. For the comparison 
the propagation of incident pulse trough the same quasi- 
harmonic wave at 72 = 73 = 0,01 is shown in Fig 10 
There is no fundamental difference between this case and 
the case of more high frequency quasi-harmonic wave Fig. 



VII. CONCLUSION 

The propagation of electromagnetic solitary wave in 
zigzag positive-negative coupled waveguide array is inves- 
tigated. Coupling between both nearest and next near- 
est neighboring waveguides is the characteristic feature of 
this waveguide system. Alternate positive and negative 
refractive indexes result in the gap in linear waves spec- 
trum. That is dissimilarity from a convenient waveguide 
array. 

The nonlinear waves were considered under assump- 
tion that NRI waveguides are linear ones and PRI waveg- 





FIG. 10: Steady state pulses crossing the quasi-harmonic plateau 
like wave, uu^jg = 0.7, 72 = 73 = 0.01. 



uides are manufactured from the cubic nonlinear dielec- 
tric. Furthermore, we took hypothesis that the phase 
difference between neighboring waveguides is constant. 
It allows to reduce the infinite system of equations to 
system of two equations and to find the steady state so- 
lutions of them. 

The characteristic coupling constant for ANOWZA can 
be controlled by the angle between the lines connecting 
the centers of neighboring waveguides i^ij. The numerical 
simulation of the collision between two steady state soli- 
tary waves shows that complicated picture of the inter- 
action is depended on coupling constants which are char- 
acterized coupling between the next nearest neighboring 
waveguides. Thus the angle '^5 is control parameter for 
nonlinear waves in ANOWZA. 
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FIG. 8: Collision of two steady state pulses in the case of 72 
73 = 0.2. 



FIG. 11: Steady state pulses crossing the quasi-harmonic plateau 
like wave, cj^^ = 8,72 = 73 = 0.01. 
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